AMARI-CHENTSOV CONNECTIONS AND THEIR GEODESICS 
ON HOMOGENEOUS SPACES OF DIFFEOMORPHISM GROUPS 



JONATAN LENELLS AND GERARD MISIOLEK 

Abstract. We study the family of a-connections of Amari-Chentsov on the 
| homogeneous space T>(M)/V^(M) of diffeomorphisms modulo volume-preserving 

diffeomorphims of a compact manifold M. We show that in some cases their 
geodesic equations yield completely integrable Hamiltonian systems. 
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1. Introduction 

O: 

In a recent paper [5] several key notions of geometric statistics were developed 
on diffeomorphism groups and their quotient spaces equipped with right-invariant 
metrics. For example, the metric defined by (the homogeneous part of) the Sobolev 
H 1 inner product of vector fields on the underlying compact manifold was shown 
to induce on the quotient of the diffeomorphism group by its subgroup of volume- 
preserving diffeomorphisms an infinite-dimensional analogue of the Fisher-Rao met- 
ric while geodesies of its Levi-Civita connection - when the underlying manifold is 
the circle - were shown to be related to solutions of a well-known one-dimensional 
completely integrable equation [6l[7]. Furthermore, the authors also described ana- 
logues of the so-called a-connections introduced in geometric statistics by Chentsov 
[3] and Amari [1] and pointed out integrability of another geodesic equation corre- 
sponding to a = — 1. 

Our goals in this paper are to provide a proof of the construction in [5] of Amari- 
Chentsov connections for circle diffeomorphisms (Theorem 12. ip . to generalize this 
construction to diffeomorphism groups of higher-dimensional manifolds (Theorem 
13. ip and finally, as a by-product, to show integrability of the geodesic equations 
corresponding to a = 1 (Theorem 12. 6[ Corollary 13. 3p . 

1.1. Diffeomorphism groups and the Fisher- Rao metric. Let M be a com- 
pact Riemannian manifold without boundary. Let T>(M) denote the group of 
smooth diffeomorphisms of M and let P M (M) be its subgroup of diffeomorphisms 
preserving the volume form fi on M. It is well-known that the completion of T>[M) 
(resp. X> M (M)) in the H s Sobolev norm with s > n/2 + 1 can be equipped with 
the structure of a smooth Hilbert manifold whose tangent space at the identity 
diffeomorphism e consists of all H s vector fields (resp. all divergence-free H s vector 
fields) on M. However, for what follows it will be sufficient to work with the smooth 
category. It will also be convenient to normalize the Riemannian volume /x(M) = 1. 
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For any rj G V(M) and any V, W G T V V(M) we set 

(1.1) (V,W)hi = \ divvdivwdfi 

4 Jm 

where V = v o rj and W = w o rj with v, w G T e T>(M) to obtain a right-invariant 
(degenerate) H 1 metric on T>{M). 

The geometry of this metric turns out to be particularly remarkable. The ho- 
mogeneous space X>(M)/X> /J (M) can be naturally identified with the set of smooth 
probability densities, i.e., smooth functions p > on M satisfying the condition 
j M pdp =1. If iV is an n-dimensional submanifold of such densities p = Pt u ..t n 
parameterized by (ii, . . . ,t n ) GK™ then recall that the Fisher-Rao metric on N is 
given by the formula 

^ f dlogpdhgp . 

(1.2) 9ij = J M - - pdp l< hJ <n. 

It turns out that the right-invariant H 1 metric defined by (11.11) on T>(M) descends 
to a (nondegenerate) metric on V(M)/V^(M) and the map rj i-> */ Jac^rj defines an 
isometry between T>(M) /V^M) and a subset of the unit sphere in L 2 (M, dp) with 
the canonical round metric. Furthermore, the restriction of the H l metric to any 
finite-dimensional submanifold iV of P(M)/T' M (M) coincides with the Fisher- Rao 
metric (11.21) on iV while its Riemannian distance is the spherical Hellinger distance 
between probability densities on M. Proofs of these statements can be found in 
Section 3 of [S]. 

Thus, in the framework of diffeomorphism groups, information geometry associ- 
ated with the Fisher-Rao metric and its spherical Hellinger distance can be viewed 
as an H 1 analogue of standard optimal transport associated with the metricQ on 
X>(M)/D M (M) induced by the (non-invariant) L 2 metric on T>(M) and whose Rie- 
mannian distance is the celebrated Kantorovich (or Wasserstein) distance, cf. [9]. 

1.2. Divergence functions and dual connections. Recall that a divergence on 
an n-dimensional manifold iV is a smooth function D : N x iV — > R satisfying 
D(p\\q) > with equality if and only if p = q and such that the matrix defined 
in a chart at p G iV by 

(1-3) 0§(p) = -^D{p\\q)\ p=q 1 < i,j < n 

is strictly positive definite for every p G N. Equation (II. 3p defines a Riemannian 
metric on N with covariant derivative determined by 

(1-4) rg fc = -^i-D(p\\ q )\ p=q 1 < i,j,k< n. 

In what follows we shall consider on T>(M) x T>(M) the functions 



(1.5) Z>< a >(eito) = (l - ^(Jac M e)^(Jac^)^^ , 

(1.6) D^ 1 \av) = D {1 \vU) = ] I (logJac^-logJac^)Jac^du 

4 Jm 



-1 < a < 1 



lr rhis metric is sometimes called Otto's metric. 
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These functions are well-defined on the homogeneous space T>(M) /T>^{M) and sat- 
isfy D^ a \^\\r]) > with equality if and only if £ and rj project onto the same 
probability density on M. 

Furthermore, recall that two affine connections V and V* on a Riemannian man- 
ifold N are called dual (or conjugate) relative to the Riemannian metric if for any 
vector fields U, V and W they satisfy W(U,V) = (V W U,V) + (U,V* W V). Basic 
facts about dual connections and divergences can be found in p] or p|]. 

2. The one-dimensional case: V(S 1 )/Rot(S 1 ) 

We first consider the case when the underlying manifold is the circle S 1 = R/Z. In 
this case D^S 1 ) is the space of rigid rotations Rot(S' 1 ) ~ S 1 . It will be convenient 
to identify the homogeneous space V(S 1 ) /Ko^S 1 ) with the subgroup of circle diffeo- 
morphisms which fix a prescribed point, for example with {rj G ©(S 11 ) : 7/(0) = 0}. 
Its tangent space at the identity can be identified with the space of smooth periodic 
functions vanishing at 0. Furthermore, given any such function u(x) we define the 
operator 

i-x ry rl ry 

(2.1) A~ l u{x) = — / / u(z) dzdy + x / / u(z) dzdy 

Jo Jo Jo Jo 

i.e., the inverse of A = —d%. 

The following is a reformulation of Theorem 6.2 stated (without proof) in [5]. 
Our first objective is to provide a proof of this result. 

Theorem 2.1 (cf. [5\, Section 6). 

(i) Each divergence induces on D(S' 1 )/Rot(S' 1 ) the H 1 metric and an affine 
connection whose Christoffel symbols are given by 

(2.2) T<f\W, V) = -l±^{A~ l d x [(V o ri-'UW o r,- l ) x ] } o v , -1 < a < 1. 

(ii) For any a the connections and V ( - _Q ' ) are dual with respect to the H 1 
metric and is the self-dual Levi-Civita connection. 

(iii) The geodesic equation o/V^ onT)(S 1 ) /Rot(S' 1 ) is the generalized Proudman- 
Johnson equation 

(2.3) u txx + (2 - a)u x u xx + uu xxx = 0. 

In particular, a = yields the completely integrable Hunter-Saxton equation 

(2.4) u txx + 2u x u xx + uu xxx = 

and a = — 1 yields the completely integrable fi-Burgers equation 

(2.5) 

^txx ~t~ 3u x u xx ~\~ wu xxx 0. 

Remark 2.2. The equation corresponding to a = 1 also turns out to be integrable 
and its solutions can be given explicitly, see Theorem 12.61 below. 

Proof. The metrics induced by and their connections can be calculated es- 
sentially as in finite dimensions using formulas ( II. 3p and ( ll.4p . We first assume 
that a 7^ ±1. Given any tangent vectors V and W at ry G V^S 1 ) let r)(s,t) be 
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a two-parameter family of diffeomorphisms in ©(S 11 ) such that 77(0, 0) = 77 with 
J^(0, 0) = V- and ' 



9 77(0, 0) = V and J^O, 0) = W. From (OJ and (HZ} we have 



^W0- = -iUSU^ (B) (^0)|K0,t)) 



( 2 -6) =T-^£U4U/ ^,0)^(0,*)^ 

s 1 



1 

1 ~a a 

1 /■ ^1+2 1 f V X W X 

= 7 / Kb VPs 77s 2 77^ 2 dx = - dx 

4 751 4 J s i r] x 

Now suppose that W is a vector field on VIS 1 ) defined in a neighborhood of r\. 
Let r)(s,t,r) be a three-parameter family of diffeomorphisms such that 7/(0,0,0) = 
77 with J^(0, 0, 0) = V, §f]{ s -> O5 0) = W^(«,o,o) f° r & U sufficiently small s, and 

7/(0, 0,0) = Z. It is clear that such a map r)(s,t,r) exists. From the formulas 
(Oil . (USD and dUBj) we have 



9s I s =0 dt \t=0dr \r=0 

' 9 1 a 1 a 



D^(r/(M,0)||r/(0,0,r)) 



1 _ q;2 3s I s=0 dt 1 t=0 3r I r=0 



7/^(5, t,0) 2 a ^(0,0,r) 2°da; 



s 1 



dx 



(2.7) =-^jL\ s=Q J ^ W x ( V (s, 0, 0))tfe(s, 0, O)^^;^^ 

= - / (DW ■ V) x Z xri - l dx - — — / W X V X 7] X 2 Z xVx 2 

= i / {« DW • r ) v' 1 )* - l -^r^ Vo o tr 1 ).}^ o v -% 

J s 

and integrating by parts and using the fact that Z is arbitrary we find that 

(yPw) v = (dw ■ v){r,) - r$(w, v) 

where the Christoffel map is given by ( 12. 2D . 

We will use the same notation for calculations in the remaining cases a = ±1. 
From ( 11 .ip and ( 11.61) we have 

iy,w)- x = (v,w)i = -f s \ s=o §- t \ t=o D^\ v (s,o)\\ v (o,t)) 

(2-8) =-^i|^SU/ i »b(«,0)(log^( aj 0)-bg^(0 l t))da ; 

4 JS 1 Vx ^ Js 1 Vx 

= (V,W) H1 . 
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The corresponding affine connections can be obtained as in (12. 7p using (11.41) . ( II. 6p 
and (I2.8p . When a = — 1 we have 

= 4L=of Lo / rix(s,t,0)-^dx 
Js 1 Vx 

= sU/ ^M,o)f^ 

= / (£W-V%— dz 
from which we deduce that 

(2.9) r^(W,V) = 0. 
When a = 1 an analogous calculation gives 

(2.10) rW(W, V) = -A- 1 ^ {(V o rT 1 )^ o r/" 1 ),)} o r/. 

This establishes the first part of the theorem. 

For the second part we need to verify that for any vector fields X, Y and Z on 
P(5 1 )/Rot(5 1 ) we have 

(2.11) X(Y, Z) H1 = <v£V, Z) H1 + (Y, V ( ~ a) Z) Hl . 

This can be done either by a direct calculation as above or else it can be deduced 
from general properties of divergences of the type (11. 5ft and ( 11. 6ft which are discussed 
in Chapter 3 of pQ. The fact that V^ -* is the Levi-Civita connection of the H 1 metric 
follows at once from ( 12. lip . 

The equation for geodesies of V < - Q - ) on V(S 1 )/Kot(S 1 ) has the form 

dt 2 7 V dt ' dt)' 

Setting dj/dt = M07 defines a time-dependent vector field u on the circle S 1 (i.e., a 
periodic function vanishing at x — 0). Differentiating this relation with respect to 
t and eliminating the first and second derivatives of 7 from the geodesic equation 
gives 

(u t + uu x ) 07 = T { ^\u 07,110 7). 

Using (I2.2p and composing both sides with 7 -1 we obtain a nonlinear pseudodiffer- 
ential equation 

u t + uu x = —^-A^dxfal) 

which we can rewrite nonlinear PDE 

yielding (EH]). □ 

Remark 2.3. The Hunter-Saxton equation (12.41) can be alternatively derived by 
observing that it is is the Euler- Arnold equation of V*- - 1 on T e V(S 1 ) /Kot(S 1 ) and 
as such it is obtained from the geodesic equation of the right-invariant H 1 metric 
(11.11) by a standard reduction procedure, see j6]. 
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Remark 2.4. Another form of the Proudman- Johnson equation can be obtained 
by integrating (12.31) in the x variable 

1 — a 2 
u tx + uu xx H — u x = C (t) 

where C{t) = — -4p J gl uldx. Observe that C(t) is a conserved integral of the 
equation (12.4p when a = 0. 

Remark 2.5 (a-curvature). Using the Christoffel symbols (12. 2p it is possible to 
calculate the curvature of the a-connections. It turns out to be proportional to the 
curvature of the H 1 metric, i.e. for any vector fields X, Y and Z on V(S 1 )/Rot(S 1 ) 
we have 

(2.12) R("\X, Y)Z = (1 - a 2 ) (x(Y, Z) til + Y(X, Z) A i) ■ 

This formula can be computed as in finite dimensions; see [8] where a different 
choice of parameters is made. 

It turns out that the geodesic equation corresponding to (12. 3 p with a = 1 can 
be integrated as well, albeit indirectly, by constructing afline coordinates for V^. 
Observe that from (" 12 . 1 2[) we already know that the connections V*- -1 -* and V*- 1 ** are 
flat. In the former case this is also evident from (12.91). 



Theorem 2.6. The geodesic equation 

(2.13) u txx + u x u xx + uu xxx = 

of is integrable. Its general solution is given by 

dri _i / \ f*e a ^ t+b Mdy 

2.14) u = — o r\ where rjlt.x) — -. r- , , -. . 

y ' dt ' J s i e a ^ t+b ^dx 

and where a,b are smooth mean-zero functions on S . 

Proof. We will construct a chart on T>{S 1 ) /Rot(S' 1 ) in which the Christoffel symbols 
of vanish identically. Consider the map 

(2.15) r] h-» (j){rj) = \ogn x - \ \ogr] x dx 

Js 1 

from P(S' 1 )/Rot(S' 1 ) to the space of smooth periodic mean-zero functions. To see 
how the Christoffel symbols transform under r] i— y fj = (j)(rj) we first compute the 
derivatives 

o nm) ^-f^ dx , Dlm v) ^Y^rv^ ix 

for any V, W E T r ,P(S 1 )/Rot(S 1 ). Next, from (J2U) and fl2~T0D we obtain 
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V X W X , f V X W X ^ _ (A- l (v x w x ) x o rj) x f (A~ 1 (v x w x ) x o rf) x ^ 



-(v x w x )or}+ / (v x w x ) o r] dx - (A 1 {y x w x ) x ) x orj+ I (A 1 (v x w x ) x ) x o rj dx 
Js 1 Js 1 

-(v x w x ) o r] + J (v x w x ) orjdx - ( - v x w x + j v x w x dx^j o 77 

+ 



DitOj; + / v x w x dx^j o t] dx = 0, 



where d = 7o r/ _1 and w = o t]^ 1 . 

We can now solve ( 12 . 13[) as follows. Since r^ 1 -* = all geodesies of in the 
afflne coordinates are the straight lines which can be written as 

t — > fj(t, x) = a(x)t + b(x) x £ S 1 

for some smooth mean-zero periodic functions a and b. Thus, given any such func- 
tions to construct a general solution u it is sufficient to (i) invert the map <ft i n (I2.15P 
to obtain the flow t — > n(t) = <fi~ l fi(t) and (ii) right-translate the velocity vector of 
rj(t) to the tangent space at the identity in V(S l ) /Rot(S' 1 )E The required formulas 
are those in f)2.14p . □ 

It is worth pointing out that the above proof manifests integrability of ( I2.13P in 
that it provides an explicit change of coordinates that linearizes the flow in the same 
spirit as the inverse scattering transform formalism. 



3. The ^-dimensional case: V(M)/V fl (M) 

We now turn to the general case when M is an n- dimensional compact Riemann- 
ian manifold without boundary and work with the coset space X>(M)/X> M (M). Our 
next result is stated in analogy with Theorem 12.11 

Theorem 3.1. 

(i) Each divergence induces on the quotient D(M)/T>p(M) the H 1 metric 
f 1 1.11) and an afflne connection given for a right-invariant vector field 
W v = w o r] and a tangent vector V = v o r\ by 

(3.1) (y ( y ] W) v = - ^A^dfddivw -v + ^-^divwdivv^ or] 

where a G [— 1, 1] and A = d5 + 5d denotes the Laplace- de Rham operator. 

(ii) For any a the connections V*-"* 1 and V ( ' _a * ) are dual with respect to the H l 
metric and is the Levi-Civita connection. 

(hi) The geodesic equation o/V^ onV(M)/T) fJi (M) is equivalent to the following 
nonlinear PDE 

(3.2) dipt + dL U dip + (1 — a)<f dip = 0, <p = divu. 



In the language of fluid dynamics the second step corresponds to going from Lagrangian to 
Eulerian coordinates. 
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Proof. Despite the fact that we are now working with cosets the computations 
involved in the proof of statements (i) and (ii) are similar to the one-dimensional 
case in Section |2j We will prove (i) for a G (— 1, 1); the proofs for a = ±1 are 
analogous and will be omitted. 

Let W be a right-invariant vector field on T>(M) defined in a neighborhood of 
r] G T>(M). Given V, Z G T fl T>(M) let rj(s,t,r) be a three-parameter family of 
diffeomorphisms such that 7/(0, 0, 0) = 77 with jrr)(0, 0, 0) = V, §fr)(s, 0, 0) = W^ujafi) 



for all sufficiently small s and ^^(0, 0, 0) = Z. Using the identity 

jy t =o Ja( VK s >*> ) = div(W,( a , ,o) °^ _1 (s,0,0)) o 77(s,0,0) Jac M 7/(,s,0,0) 
= div u; o r)(s, 0, 0) Jac Ai ?7(s, 0, 0), 

where W v = w o 77 and similar computations for the partial derivatives in s and r 
we have 



(V { y ) W ) Z) a = - [ dw({V i y ) W)0 V ' 1 )dw{Z0 V - 1 )dfl 

4 J M 

r=0 D^( V (s,t,0)U0,0,r)) 



d_\ d_ I _9_ 
<9s l s =0dt lt=0 9r \r- 

1 a I 9 



(Jac /i 77(s,t, 0)) 2" (Jac^O, 0, r)) 2°^ 



7^L=o / divwor ?( s > >0) (Jac M ?7(s,0, O))^" div(Z o 77 x ) o 77 (Jac^) '^d/i 
4 Jm 

7 / i rfdi 



div w o 77 • V 



1 — 1 

H — div wot] div(V o rj~ l ) o 77 1 div(Z o r/ _1 ) o 77 Jac M 7/ d// 

- / /ddivit; • v H div w div 7j | div(Z o i]^ 1 ) d/x, 

4 y M I 2 J 



where V" = v o 77. Using integration by parts 



f div XdLt = - f5X"dLi = - df(X)dfi 

M JM JM 



and the fact that Z is arbitrary, we find 



ddiv{(V ( y ) W)o V - 1 } =df i.e. dw{(V { y ) W)o V - 1 } = f - [ fdLl, 

Jm 

where / = d div w ■ v + divwdivw. The expression in (13 .ip now follows from 
the next lemma. 

Lemma 3.2. Let A be the Laplace- de Rham operator. Then 

div(-(A- 1 d/)») = / - / fdLi, f e C°°(M). 

Proof. If g = div(— (A _1 <i/)^) = 5A~ l df then Ag = 5df = Af and since the kernel 
of A acting on functions consists of the constants, we deduce that g — f + c for 
some constant c. Integrating over M yields = J M fdii + c, which determines c. □ 
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Calculation of for an arbitrary (not necessarily right-invariant) vector field 

W can be reduced to the right-invariant case as follows. Let r) G T>(M) and let W R 
be the right-invariant vector field with W R = W n . Then 

(3.3) 

(vPw) v = (vP(w - w R )) v + (vPw R ) v = (c v (w - w R )\ + (v{?V fl )„, 

where Ly denotes the Lie derivative in the direction of any vector field V such that 
V v = V. Indeed, if Z is a vector fields satisfying Z v = 0, then locally 

(V V Z) V = DZirj) ■ V - T V (Z V , V) = DZ( V ) ■ V = [V, Z} v = (C V Z) V . 

Using ( 13. ip and (13. 3p we can compute the affine connection for arbitrary fields 
on T>(M) and this completes the proof of (i). 

As in the proof of Theorem I2.1[ (ii) can be established by a direct calculation or, 
alternatively, it can be deduced from the general properties of divergences of the 
type (OD and (Oil . 

Regarding (iii) let j(t) be a curve in T>(M) with 7(0) = e and 7 = dj/dt = u o 7. 
In the appendix we prove that (13.31) implies 

(3.4) V^7 = Mt o7 + vfy. 

Using (13.11) and (13.41) the geodesic equation Vjy 7 = can be written as 

1 — a 

Ut = (A~ 1 dfY where f = ddivu-u-\ — (div-zz) 2 

so that setting (p = div u and taking the exterior derivative we get 

—dipt + Sdu t = di u dip + (1 — a)(pd<p. 

The relation u\ = A~ 1 df implies Adu\ = dAu\ = hence, in particular, 5du\ = 0. 
This proves ( 13. 21) . □ 

Finally, we turn to integrability of the geodesic equations in ( 13 .21) which to the best 
of our knowledge has not been studied in the literature before, with the exception 
of the case a = in [5] . 

Corollary 3.3. The geodesic equations o/V^ on T>(M) /T)^(M) corresponding to 
( 13. 21) with a = 0, ±1 are integrable equations in any dimension n. 

Proof. Equation (13. 2p with a = was derived and shown to be integrable in [5]. 
Since and V*- -1 ^ are flat, integrability of the other two equations can be estab- 
lished similarly to their one- dimensional analogues in Section [2j We will consider 
the more complicated case a = 1 with the corresponding equation 

(3.5) dipt + di u dip = 0, ip — divu. 
Proceeding as in the proof of Theorem 12.61 we first show that the map 

rj 1— > (f)(rj) = log Jac^r/ — / log Jac^ d^i 

Jm 

from T>(M) /T>^(M) to the space of smooth mean-zero functions on M defines an 
affine chart for V^. Thus, as before, the geodesies of V 1 - 1 ^ in the affine coordinates 
defined by are the straight lines 

t — y fj(t, x) — a(x)t + b(x) xeM 
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where a and b are smooth mean-zero functions on M. From this we find 



3ac u ri(t,x) = -7 ,,,,,,,, x E M 



and combining this expression with the identity J^Jac^r/) = (^077) Jac^r/, we obtain 

f..n( 

(3.6) (p(t,rj(t,x)) = a(x) - 



J M a(x)e a( ^ t+b(x Un 



j M ea(x)t+b(x) dfi 

Note that the time derivative of <p o rj is independent of x, so that 

<kp t + dt u d(p = d((<f o T)) t o r/" 1 ) = 
which shows that (13. 6p solves the equations in (13. 5p . □ 

Appendix A. Proof of ( 13.41) 
Let 7 (t) be a curve in V(M) with j(t) — u(t) o j(t). For any t > we compute 

Ut(*o) l(to) = I ^ (7(0 7(0~ X ) } ° 7(*o) 
(A.l) =7(*o)+£>(7(*o))-{| 

But differentiating the relation 7(t) _1 o 7 (t) = id with respect to i we find 

lit)' 1 } o 7 (t ) + J D(7(to)" 1 ) ° 7(*o) • 7(*o) = 0. 



7(*)~ 1 |> °7(*o)- 

t=t 



t=t 



d 
dl 

Thus, (TA~Tj) yields 

«t(*o) o 7(*o) = 7(*o) - £(t(*o)) ■ D^ito)- 1 ) o 7 (t ) ■ 7 (t 
d 



{ 7 (t) - 7 (t ) o 7 (t )™ 1 o 7 (t)} = D( 7 - 7 R ) ■ 7(*o)- 
t=t 



Equation (13. 4p now follows from (13. 3p . 
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